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Birkhoff's theorem is discussed in the frame of f(R) gravity by using its scalar-tensor represen- 
tation. Modified gravity has become very popular at recent times as it is able to reproduce the 
unification of inflation and late-time acceleration with no need of a dark energy component or an 
inflaton field. Here, another aspect of modified f(R) gravity is studied, specifically the range of 
validity of Birkhoff's theorem, compared with another alternative to General Relativity, the well 
known Brans-Dicke theory. As a novelty, here both theories are studied by using a conformal trans- 
£_) ' formation and writing the actions in the Einstein frame, where spherically symmetric solutions are 

studied by using perturbation techniques. The differences between both theories are analyzed as 

^*\ ■ well as the validity of the theorem within the Jordan and Einstein frames, where interesting results 

are obtained. 
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In the last decade, modified theories of gravity have become very popular as they can explain the accelerating 
_». | expansion of the Universe without need to introduce dark energy. Such theories can even provide an alternative 
^. ■ description of inflation without any need for an inflaton field. Hence, dark energy and inflationary eras appear as 
a natural consequence of the law of gravity, which has a different behavior depending on the scale. Then, gravity 
is supposed to be well described by General Relativity at local scales, while the additional terms introduced in 
theories of modified gravity become important at cosmological scales, so that the cosmological evolution can be 
described in purely gravitational terms with no need to assume additional components (a seminal paper in this sense 
is |l|). Special interest in this gravitational paradigm is related with f{R) gravity (for general reviews, see [3, |3|), 
which has a simple structure, where more complex functions of the Ricci scalar than in the Hilbert-Einstein action 
are considered, but with no other curvature invariant or non-local ternr. Within f(R) gravity, reconstruction of 
inflationary and dark energy epochs can be easily performed and useful techniques have been developed in order 
to reconstruct the appropriate action that can describe the Universe evolution (see [4]). Also, the so-called viable 
models of f(R) gravity can avoid violations of local gravity tests and matter instabilities (see Ref. 5]). In this sense, 
the study of spherically symmetric solutions in the frame of f(R) gravity turns out an important and interesting 
question, which has been already dealt in Ref. |6| , as well as the study of the Newtonian limit of the theory, which 
can constrain the form of the action (see Ref. 7]). In most cases the reconstruction of the action f(R), specially when 
dealing with Friedmann-Lemaitre-Robertson- Walker metrics, is done by using an auxiliary scalar field, since f(R) 
gravity is equivalent to a kind of Brans-Dicke theory with a non-propagating scalar field and a non-null potential. 
However, both theories seem to exhibit a different behavior in the study of perturbations, as it was pointed out in 
Ref. |8}. This is probably due to the high non-linearity exhibited in f(R) gravities, which can be alleviated when one 
deals with its equivalent picture in a Brans-Dicke-like theory. 

At the present paper, we are interested to study spherically symmetric solutions and specifically Birkhoff's theorem 
in f(R) gravity and its scalar-tensor representation. It is well known that, in f(R) gravity as in Brans-Dicke theory, 
Birkhoff's theorem is not satisfied unless strong restrictions are imposed on the scalar curvature and on the scalar 
field respectively, Ref. [9]. However, here we are interested to study the range of validity of Birkhoff's theorem by 
performing perturbations around a background solution in f(R) gravity, where an auxiliary scalar field is used, and 
compare the results with those for Brans-Dicke theory. In order to obtain a real comparative of the perturbations 
in both theories, we shall use the mathematically equivalent picture defined in the Einstein frame, which is related 
with the original one, usually called Jordan frame, by a conformal transformation. Then, by writing /(i?) gravity 
and Brans-Dicke theory in the Einstein frame, we study the perturbations around a background solution, which can 
be seen as a new technique to compare both theories. The results obtained in the Einstein and the Jordan frames 
are studied, where it is found different information on the Birkhoff's theorem for each frame, what may suggest the 



non-physically equivalence between both, as it was already pointed out in Ref. [10( in a cosmological context. 

Hence, the paper is organized as follows: next section is aimed to the framework of the paper, where f(R) gravity and 
its scalar-tensor representation (the so-called O'Halon theory) and Brans-Dicke theory are introduced. The relation 
between the Jordan frame, where both theories are defined, and the Einstein frame is briefly introduced in Section 
III, where the transformation of a general spherically symmetric solution is obtained. In section IV, the perturbative 
study of Birkhoff 's theorem in the Einstein frame is done, where the main results of the paper arc obtained. Finally 
in section V, we summarize and present some conclusions. 

THE FRAMEWORK 

Let us start writing the action and field equations that describe a general f(R) theory. We consider the action, 

S = Jd 4 x^[f(R) + 2 K 2 £ m ] . (1) 

where k 2 = 8ttG, R is the Ricci scalar and £ m stands for the Lagrangian corresponding to matter of some kind. Note 
that this action reduces to General Relativity for f(R) = R, and the corresponding field equations turn out of second 
order, while for the more general action ([1]), the equations are fourth order, which are obtained by varying the action 
(UJ respect to the metric tensor g^, 
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Here the primes denote derivatives with respect R, and the energy-momentum tensor is given by X)}„ = 
By taking the trace of the equation (J5J , one obtains an extra equation given by, 

3af(R) + f'(R)R-2f 
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Due to the extra terms in the field equations, these theories can reproduce quite well the cosmic acceleration with 
no need to introduce a dark energy component or a cosmological constant, and even the inflationary epoch can be 
described as a pure gravitational effect (for some literature, see Refs. [4) and 5]). The study of this class of theories 
in contexts apart from cosmology is a fundamental issue to understand these theories and to provide new results. At 
the current paper we are interested to study spherically symmetric solutions (see Ref. [6J]-[8|]) and the range of validity 
of the Birkhoff 's theorem in f(R) gravity. We will explore the relation between f(R) gravity and scalar-tensor theory 
in the Jordan and Einstein frames, and we compare the results obtained for each frame. For that, let us write the 
general action for a Brans-Dicke-like theory, 



Sbd = d xy^g 



cbR - -g^V^,4> ~ V(4>) + 2k 2 C T: 



(4) 



Here we assume w to be a constant. This action is written in the so-called Jordan frame, which is related with the 
Einstein frame by means of a conformal transformation. The field equations for this action are obtained by varying 
the action respect the metric tensor g^ and the scalar field 4>, 
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d<p 
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The action (UJ) becomes the Lagrangian describing O'Halon theory for a null kinetic term w = 0, such that the actions 
([1]) and (|4]) turn out equivalent. This is straightforward to see, just by varying the action ((4]) with respect the scalar 
field 4>, where the relations f'(R) = <j> and V"'(</>) = R are obtained, which reflects the extra scalar degree of freedom 
that f(R) gravity owns, where the trace equation ([3]) can be seen as a Klein-Gordon equation. On the other hand, the 
original Brans-Dicke theory is given by taking a null scalar potential V{<j>) = and w ^ in the action (|4]). Despite 
the similarity between the two theories, when the Newtonian limit is performed, the well known PPN parameter 



for Brans-Dicke theory, given by jbd = ^+^jj does not correspond to the PPN parameter for the O'Halon theory 
(equivalent to f(R) gravity) by taking w = 0, which gives 7 = 1/2, but the PPN parameter for f(R) gravity is given 
by a more complex expression as it is shown in Ref. [8|. Nevertheless, here we are interested to study perturbations 
around a given background solution, particularly we will study spherically symmetric solutions of the type, 

ds 2 = -A(r, t)dt 2 + B(r, t)dr 2 + r 2 dVL , (6) 

where dQ. is the metric of a 2-sphere. The metric ^ can be considered as the most general spherically symmetric 
solution, where the coordinates have been chosen to avoid crossed terms between the spatial and time terms. For the 
theories described by the field equations (|5|), the only way to ensure a static solution ([6]), that is A(r,t) — A(r) and 
B(r, t) — B{r), is to impose a time independent scalar field <f>(r, t) = <j){r), which implies in the equivalent f(R) gravity 
representation to assume a time independent Ricci scalar in the field equations (jSJl (see Ref. [9]). At the present paper, 
we do not impose any condition on the scalar field to study spherically symmetric solutions ©, but we assume that 
both theories introduce small corrections to General Relativity in order to avoid violations of local gravity tests, as it 
has been pointed out in several works for f(R) gravity (see Ref. [5|). We are interested to study perturbations around 
a background solution in Brans-Dicke theory and f(R) gravity (by means of its O'Halon description) via the Einstein 
frame, where the action (j4|) can be transformed by performing a conformal transformation as it is well known. This 
also serves to compare the results within the Jordan and Einstein frames, where the possible non-physical equivalence 
between both frames is shown. 

RELATION BETWEEN f(R) GRAVITY AND MINIMALLY COUPLED SCALAR-TENSOR THEORY 

VIA CONFORMAL TRANSFORMATIONS 

As it was pointed out in Ref. [8J], the presence of a self- interacting scalar potential in O'Halon theory compared with 
BD theory, makes the perturbations in the Newtonian limit of a spherically symmetric solution to not be comparable, 
so that here we use the mathematically equivalent Einstein frame to study both theories and compare the results 
between the different frames. 

In the Einstein frame, the scalar field in the action is minimally coupled to the gravitational field, and the equations 
acquire a shape similar to that of General Relativity. The relation between the Jordan and Einstein frames is given 
by the conformal transformation, 

gE lw = ti 2 gnv, where n 2 = . (7) 



Then, the action (J3J) in the Einstein frame becomes 

2w + 3„ V{4>) 



S E = d Xy/-gE 



R * -^S-V^V"* - -W- + 2k2£ i 
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where the subscript e denotes that the variables are defined in the Einstein frame and the Lagrangian of matter is 
given by £^ m = iAn ( -zsQeimv ) ■ Note that O'Halon theory is described in the Einstein frame by the action © 
with w = 0, while Brans-Dicke theory is defined by V((f>) = and w 7^ 0. In order to simplify the equations, we can 
redefine the scalar field as (j> — e^/^ 3 " 1 " 2 ™, and the action (j8]) takes the form, 

1 



Se = d 4 xy/-g E 



Re - ^V^V> - U(<p) + 2 K 2 d 
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here the scalar potential is defined as U(ip) — e v ^ 3+2w V ((j)(ip)) . The field equations can be obtained by varying the 
action ([9]) with respect to <7_e m „ and <fi, 

Re^v - ~9eij. v Re = -xd^ipducp - -gEp,u [d a <pd a cp + U(<p)} + K 2 T%£l , (10) 

Dv -*M = -2K* s(c s Bm) , (ii) 

dip dip 
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where the energy-momentum tensor is defined as T E ' = ._ °rf g 1 . By taking the trace of the first equation in 
(|10p . the auxiliary equation is obtained, 

R E = \d<,yd°y + 2U{tp) + n 2 T% n) . (12) 



Hence, we will use the equations (ITU1) - (|12I) to study the perturbations created by the scalar field ip in the Einstein frame 
around a background solution. Nevertheless, the spherically symmetric solution ([5]) is transformed in the Einstein 
frame by the conformal factor (JT)), and it yields, 

ds% = II 2 (r, t)ds 2 = -Il 2 A(r, t)dt 2 + U 2 B(r, t)dr 2 + U 2 r 2 dQ . (13) 

By redefining the coordinates, the metric can be rewritten in the more familiar form given in ^ , taking p 2 = II 2 (r, t)r 2 , 
and t' — T(t, r) to avoid cross terms, the metric (fT3")l is rewritten as, 

ds% = -C(p, t')dt' 2 + D{p, t')dp 2 + p 2 dil . (14) 

Hence, we will study the perturbations in the metric (J14I) around a background solution, which will provide the range 
of validity of Birkhoff 's theorem in the Einstein frame in terms of the order of the perturbations, and for comparing 
the results obtained for O'Halon and Brans-Dicke theory, as well as the study the differences between the Einstein 
and Jordan frames. 

SPHERICALLY SYMMETRIC SOLUTIONS AND BIRKHOFF'S THEOREM 

Let us now study the metric (IT4"|) for the general scalar-tensor theory defined by the action ©, and explore the 
range where such metric is static. It is well known that for Einstein's field equations, the only solution in vacuum for a 
spherically symmetric metric is given by the Schwarzschild solution, or Schwarzschild-(A)dS solution if a cosmological 
constant is included in the field equations. This result, called Birkhoff's theorem, was proved independently by G. 



D. Birkhoff [ll| and J. T Jebsen J12J, and it states the following, 



Birkhoff's theorem: a spherically symmetric solution of the vacuum Einstein equations is always static in a 
region where the time coordinate remains time-like and spatial coordinates stay space-like. 

Here we want to explore the range of validity of this theorem for f(R) gravity and Brans-Dicke theory trough out 
their description in the Einstein frame given by the action © and the field equations (fj"0|) - ([i"2j) . It was pointed out 
in Ref. [7| that in f(R) gravity, the Birkhoff's theorem is valid up to fourth order in perturbations in the Newtonian 
limit of the theory, what implies the weak field regime. Nevertheless, for other regimes, where the scalar curvature 
is sufficiently large, such limit would not be valid, and the solution can change. Here, we consider a more general 
context, where a background solution for the metric (|14[) and the scalar field is considered, and perturbations around 
the zero-order solution are studied. Then, the metric and scalar field can be written as, 

9E, v = 9 Z + 9 { Zu, ^ (0) +^ (1) . (15) 

Here g E refers to the zero-order solution, while the perturbations are represented by g Euv - By taking the metric 
(fT4")l . the components of the metric can be written as 

gm>t> = -C(p,t>) ~ -C*°>(p,tf) - CU(p,t>) 
g Epp = D(p,t>) ~ D(°>(p,tf) +DW(p,t>) 

9E68 = P 2 

V 9E4>4, = P 2 sin 2 9 



< 



while the inverse metric, defined as usually by g E ° ' gEa\ = <5", yields, 

n t't> _ i i__ , C^\p,t') 

Ve C{p,t') — CW(p,t>) ^ (C<°)(p,t')) 2 

n PP _ 1 1 D {1) (p,t') 

9e — D(p,f) — D(<»(p,t>) (_D<°)( p .t'))2 (17) 

9e - -p* 
n iM> — l 

K ME — p* sin^ 8 

The Ricci tensor can be written in its standard form as, 

Re = g E R-Epv = g E ( r M „ iQ — ^ pa ,u + r '^^r '^ — T^ p T ua \ , (18) 



where the Christoffel symbols are defined by, 

TJL = l^B (9Ena,u + S£W,m - 9Efj,iy,cr) ■ (19) 

By the perturbed metric defined in (|16[) and (1171) . we can split the Christoffel symbols into two orders of perturbations, 
T" = r(°) + rW , as well as the Ricci scalar Re = R E + -4 ■ The perturbations also act on the scalar potential 
U(tp), such that the potential can be expanded around a background solution for the scalar field ipo, 



£%) = ^^U>-V> (0) - (20) 



Jjn{0) 
n 

By inserting the above expressions in the field equations (|10l) - (|12[) . we can split the equations into the different orders 
of perturbations. We are interested in vacuum solutions, Tg™ v — 0, and we assume that the background solution is 

given by a constant scalar field tf(°>(p, t') = ipo- The reason for this is that it is assumed that at zero-order the results 
for General Relativity has to be recovered, and the additional terms in the action are involved at higher orders only. 
Then, the equations (TIT))) and (11 1|) at zero-order are given by, 

4t-^4 0) +.4°>.o (2i) 



dif> 



, (22) 



where the cosmological constant is defined as A = \U$. The equations (|2ip are exactly the same as the Einstein field 
equations with a cosmological constant, and the solution is the well known Schwarzschild-(A)dS metric, that gives 
the zero-order solution, 

C^{p) = [D^{p)]^ = l- 2 -^-±p\ (23) 

p 3 

where p is an integration constant. Then, at zero-order we have a static solution which satisfies Birkhoff's theorem 
as expected. At first linear order, the equations (|10 p - (fT2"|) are, 



D^ 1 ) = C/^ (1) , (25) 

RW = 2£/y ^ . (26) 

Introducing the results obtained at zero order (|2"2"1) . where we have U' Q — 0, and R^ — 2Uq, the equations (|24p ~ ([2^1) 
turn out much simpler, where the field equations for the metric (|24p yield, 

4V^off£i,=0. (27) 

Note that the scalar field decouples from the metric at linear order, which is not expected at higher orders. The 
expression (|27[) is a linear system of differential equations, where the coefficients are given in terms of the zero order 
(|23p . and whose solution exhibits a very complex expression (for a complete analysis on the perturbations of a black 
hole in f(R) gravity, see Ref. [13j). Nevertheless, we can study the system in order to obtain the enough information 
about the form of the metric components at first linear order and their staticity. From the ri-equation, 

^ = ^=0 -> g M p = D U( P ,t>) = DU(p). (28) 

Hence, the rr-component of the metric is time independent. By deriving the #0-equation respect the time, it gives, 



It is straightforward to show by the rr- and tt-equations that a(t') remains an undetermined function of the time 
t', so that it can be taken to be constant a(t') = a, while the functions D^'(p) and x(p) are solutions of the 
system of differential equations (l27l) . In general, one can not obtain a general solution for D^(p) and x(p), but 
numerical methods are required. However, the solutions (|28H29p provide the enough information to analyze the 
metric components. As they are time-independent, Birkhoff's theorem is satisfied at this order by the solution, and 
the metric remains static in vacuum for a general scalar-tensor theory described by the action © in the Einstein 
frame. 

At this point we have to remark the differences between O'Halon and Brans-Dicke theory described both in the 
Einstein frame by the action ©. It is straightforward to show that for BD theory, where the scalar potential 
U((p) = 0, the solutions become quite different compared with O'Halon theory. At zero-order, the field equations (f2"Tj) 
are different for each theory, while in BD the cosmological constant vanishes and the zero-order solution (|23[) turns 
out the classical Schwarzschild solution, for O'Halon theory the metric is given by Schwarzschild-(A)dS solution at 
zero order. At first linear order, the equations (pMf reduces to R Eilv — in BD theory, whose solution is also static, 
but naturally it gives a different metric compared with O'Halon theory, as it was expected. Note that up to this 
point the only term used to compare both theories is the scalar potential, while the parameter w does not play any 
role in the Einstein frame, but on the inverse conformal transformation to turn back to the Jordan frame. Then, this 
provides a new method to compare both theories in perturbations using the Einstein frame representation. In this 
case, we have found that Birkhoff's theorem is satisfied in both approaches at least to first linear order, although the 
solutions for the metric arc different for each approach. 

Let us now transform the metric components obtained in the Einstein frame (|28H29p to the Jordan one, and check 
how the metric is affected. In order to obtain the conformal transformation 0, we have to solve the scalar field 
equation at first linear order (1251) . which can be extended, 

This equation clearly has a solution where we can split up the time and radial parts of the scalar field, 

^(p,0 = ^ 1 VS 1) . (31) 

And the equation ([50)1 takes the form, 

C<°Vf) - (c*» + H<7«») 4fU _ DJJVW _ ^4l! - ■ (32) 



Here, the primes denote derivatives with respect to the radial coordinate while the dots denote with respect to the 
time coordinate. Hence, the corresponding time solution of the equation has to hold, 

% = k -+ <p$> = de™ + C 2 e-^ kt ' , (33) 

where C\2 an d k are constants. Note that the sign of k may play an important role on the stability of the background 
solution in the Jordan frame, as an exponential function is obtained for a positive value of k, while an oscillator is 
obtained for a negative k. The radial part of the solution ip p is obtained by solving the differential equation, 



c( o y ,(1) (p) 

Hence, the conformal transformation yields 



C"(o) + H(7(0) 
r 



^ 1 Hp)-{^ + U^^ 1 Hp)=0. (34) 



n 2 = 4> = ^ + 0« = &/>&&* ~ ^°) (l + l , r w) , (35) 

V v3 + 2w / 

where 4>^ — e v /V3+2™ _ Recall that O'Halon theory is defined by w — 0, so the inverse of the conformal transfor- 
mation ([35]) in BD and O'Halon theories differs only in a numerical factor, and the same differences between both 



theories that were found in the Einstein frame, remain in the Jordan one. However, the important result here refers to 



the differences between the Einstein and Jordan frames. By the transformation (1551) . it is clear that the perturbations 
on the metric in the Jordan frame will not be static at first linear order, and the Birkhoff 's theorem will not be valid in 
this frame, which differs completely from the Einstein one. Note also that, in general the result given in (|33|) . suggests 
that the zero-order solution will be unstable in the Jordan frame due to the perturbations introduced by the freedom 
of the scalar field. This fact suggests the non-physical equivalence between both frames, as it has already pointed out 
in Ref. [lOj by the analysis of both frames in a cosmological context. The transformed metric in the Jordan frame is 
given by, 

ds * = -^M-dt* + #?W + -/-<m . (36) 

At zero order, where 4>^ in (|3"5|) is a constant, the result in the Jordan frame gives also a Schwarzschild-(A)dS metric, 
and the Birkhoff's theorem is satisfied. Nevertheless, at the first linear order the metric is clearly not static, which 
contradicts the result obtained in the Einstein frame, and points to the different physical meaning of both frames. 
By redefining the radial coordinate r 2 = tfi(p,t')p 2 in (|36j) . and performing some variable transformations, the metric 
(|36p can be written in the more common form given in ((6|), where the temporal and radial components of the metric 
will depend on time, introducing a strong instability on the background solution. 

DISCUSSIONS 

In summary, we have studied the range of validity of Birkhoff's theorem in f(R) gravity through its O'Halon 
representation, and in Brans-Dicke theory, both expressed in the Einstein frame. The comparison with the results 
in the Jordan frame has been also analyzed. Here we have presented a method to compare the two theories in the 
Einstein frame instead the Jordan one, where both theories have to be studied separately in the perturbations in order 
to obtain consistent results, as pointed out in Ref. Q- We have seen that assuming a background solution described 
by a constant scalar field (as a limit to General Relativity), the zero-order solution gives naturally a static metric, 
the (A)dS-Schwarzschild solution. However, also the first linear order provides a metric that is time-independent 
in the Einstein frame, so that the Birkhoff's theorem is satisfied. Nevertheless, the fact that Brans-Dicke theory 
has a null scalar potential makes the solution different compared with O'Halon theory, although the theorem is also 
satisfied, at least for first linear order in perturbations. When non-linear terms are included, it is expected that the 
metric will not remain static in both theories. 

In parallel, we found that the result obtained in the Einstein frame about the range of validity of the Birkhoff's 
theorem is affected when a conformal transformation is applied, and the action is transformed to the Jordan frame, 
where the metric is not static at first linear order in perturbations. Hence, this can mean that the Jordan and 
the Einstein frames are not physically equivalents, at least for a perturbative analysis. While in the Einstein frame, 
Birkhoff's theorem is valid at first linear order, where the scalar field remains decoupled, the conformal transformation 
makes the metric depends on time when it is transformed to the Jordan frame, where the coupling with the scalar 
field is recovered, violating the Birkhoff's theorem. In addition, the presence of a time-dependence in the metric can 
affect the stability of the background solution, which was found to be stable in the Einstein frame. Then, as pointed 
out in Ref. |l0| . both frames may be mathematically but non-physically equivalent. 
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